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Topological quantization by controlled paths: application to Cooper pairs pumps 
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When physical systems are tunable by three classical parameters, level degeneracies may occur 
at isolated points in parameter space. A topological singularity in the phase of the degenerate 
eigenvectors exists at these points. When a path encloses such point, the accumulated geometrical 
phase is sensitive to its presence. Furthermore, surfaces in parameter space enclosing such point can 
be used to characterize the eigenvector singularities through their Chern indices, which are integers. 
They can be used to quantize a physical quantity of interest. This quantity changes continuously 
during an adiabatic evolution along a path in parameter space. Quantization requires to turn this 
path into a surface with a well defined Chern index. We analyze the conditions necessary to a 
Topological Quantization by Controlled Paths. It is applied to Cooper pair pumps. For more general 
problems, a set of four criteria are proposed to check if topological quantization is possible. 

PACS numbers: 85.25.Cp, 03.65.Vf, 74.50.+r, 74.78.Na 
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I. INTRODUCTION 

Using nanometer size Josephson junctions, a huge va- 
riety of superconducting quantum circuits can be made. 
These circuits are described by simple Hamiltonians in- 
volving a discrete set of quantum variables. They are 
typically the excess number of Cooper pairs hj on super- 
conducting elements and their canonical conjugate vari- 
ables Qj {[fij, Qj\ = i) are related to the quantum phases 
of the superconducting order parameters of the circuit is- 
lands. In addition, most circuits have tunable elements: 
they are control voltages on gates, or using magnetic 
fluxes, quenched quantum phases or Josephson couplings. 
These circuits are most often used to implement quan- 
tum logic^'^, where the quantum gates are controlled with 
voltage or resonant microwave pulses on some of gates or 
other tunable elements. 

The tunable elements of quantum circuits can also be 
used to generate adiabatic evolutions of the Hamiltonian 
as a function of the parameters. More precisely, the N 
induced gate charges Ugi and control phases ipi define 
a vector R = {rigi . . . (^i . . .} in the parameter space P 
of dimension N . Let i?c((R) . . . and |q:(R)) ... be the 
eigenenergies (bands) and eigenvectors of the Hamilto- 
nian 7i(R). At a point R and for a non-degenerate band 
a, one can construct a fiber defined by the set of all vec- 
tors {|a(R))} which may differ by a complex factor. The 
set of all these fibers defines the fiber bundle over the 
parameter space P. When the topology of this bundle 
becomes non-trivial, physical phenomena of great inter- 
est can occur. 

Parallel transport, holonomy and homotopy are cen- 
tral concepts for the physics of geometric phases. Berry's 
phase is one of these and is a relevant quantity when adi- 
abaticity conditions holds for a globally non-degenerate 
band. In this case. Berry's phase is the geometric part 
of the phase acquired by the wavefunction along an adia- 
batic cycle over a closed path in the parameter space 
P. 

Depending on the nature of the quantum system stud- 



ied, the physical consequences of the non-trivial topol- 
ogy of the eigenvector bundle are different. A number 
of physical examples have been studied in several areas 
of physics. In molecular systems, the electronic structure 
depends on the semi-classical nuclear coordinates (within 
the Born-Oppenheimer approximation) which define the 
parameter space. Their energy manifolds can have coni- 
cal intersections at isolated values of the nuclear coordi- 
nates'^"^: these so-called "diabolical points" are directly 
responsible for the change of multiplicity of rotation- 
vibrations levels as a function of nuclear coordinates. In 
molecular magnets, the magnetic energy levels depend 
on the direction and magnitude of the applied magnetic 
field (the parameter space) with respect to the molecu- 
lar axes. For some molecules, isolated degeneracies have 
also been found for specific direction and values of the 
magnetic field^'^. At these points, quantum tunnelling is 
quenched as a result of interferences caused by the wave- 
function phase changes around the "defect" . Indeed, this 
phase change takes a particular value of tt for physically 
relevant paths encircling the degeneracy, leading to the 
destructive interferences observed. 

In the few examples above, degeneracies occur at iso- 
lated conical intersections between two energy bands in 
a three-dimensional parameter space. These diaboli- 
cal points are singularities of the quantum phase field 
over the parameter space and responsible for the "ex- 
otic topologies". Closed paths, through Berry's phase, 
are sensitive probes of the topology. Closed surfaces in 
a three-dimensional parameter space are also sensitive 
to the presence of conical points through a topological 
invariant called the Chern index ci, which is an inte- 
ger number. Some of the best known phenomena in 
condensed matter physics are well-understood in terms 
of Chern indices, such as integer quantum Hall effect 
(IQHE)*'9, Thouless Pumping (TP)i'' or AC Josephson 
effect (ACJE)^°'^^. Physical quantities which can be ex- 
pressed in term of Chern indices are subject to topological 
quantization. This is why they are used in metrology. For 
instance, the IQHE gives a conversion from voltage unit 
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(Volt) to current unit (Ampere) through the resistance 
quantum = V = ; in the same way, the 

AC Josephson effect gives a conversion from voltage to 
frequency (in Hertz) through the magnetic flux quantum 

Such physics can be encountered in the simplest super- 
conducting circuit depending on three tunable parame- 
ters: the Cooper Pair Pump (CPP), where degeneracies 
occur at isolated points in the parameter space. Here, the 
parameter space is constructed from two gate voltages Vgi 
and a quenched quantum phase (p. However, an essential 
difference subsists between the examples given above and 
our problem. In the case of IQHE, TP or ACJE, the rele- 
vant physical quantity measured is directly proportional 
to the Chcrn index of the surface broTight into play (the 
magnetic Brillouin zone for the IQHE). In a CPP, one 
can only make (adiabatic) paths in the parameter space 
by modifying the parameters in order to tune the current 
delivered by the CPP. Thus, Berry's phase seems to be, 
a priori, the relevant topological quantity characterizing 
the paths. This will be shown to always be the case, 
by relating the charge transferred by the CPP through 
a path is always expressible in term of a Berry's phase, 
even for open paths. Nevertheless, Chern indices can 
also specify the value of the current for specific paths 
covering densely a closed surface enclosing a diabolical 
point. The degree of quantization of the current deliv- 
ered by this method improves exponentially with the de- 
gree of the surface coverage. This form of quantization 
is referred to as "Topological Quantization by Controlled 
Paths" (TQCP). The current will be shown to be equal 
to 2e u where v is a, characteristic frequency of the adia- 
batic cycles of pumping. For metrology, the Cooper pair 
pumping through TQCP gives a conversion between cur- 
rent and frequency: it is an effect which could be used 
to close the metrological triangle between the units of 
voltage V , current I and frequency v. 

In this paper, we emphasize what is new and specific 
to quantum circuits, and the example of the Cooper pair 
pump is an excellent case study for the concept of TQCP 
around which much of the paper is built. In Sec. II, 
we recall the topological properties of three dimensional 
parameter spaces stressing the notions of Berry's phase 
and Chern indices in the presence of degeneracies. In 
Sec. HI, the TQCP is introduced, and the computation 
method of the quantized physical quantity is given. In 
this section, the necessary criteria for TQCP are derived, 
namely 

1. The energy spectrum must be discrete (adiabatic- 
ity). 

2. The Hamiltonian depends on three continuous tun- 
able parameters, which specify the parameter space 
P. Isolated conical degeneracies between the two 
lowest eigenstates |±) must occur in P. 

3. The relevant physical observable Q (the quantity 
measured) follows Hamilton semiclassical equation 



of motion Q = {d,/H) where the parameter is pe- 
riodic. The contribution of this quantity along ge- 
ometrical paths is set by the topology of the eigen- 
vector bundle. 

4. The dynamical contributions to Q must also be 
taken into account. The topological quantization 
can be implemented only when they can be elimi- 
nated. This is possible when the (/^-dependence of 
the Hamiltonian eigenvalues can be integrated out 
using its periodicity or other symmetries of the sys- 
tem. 

In the conclusion, a full discussion of these four criteria 
is presented in light of this work. Section V is devoted 
to a practical implementation of the TQCP for Cooper 
pair pumps. Section VI shows how microwave fields can 
be used to expand the parameter space to higher dimen- 
sions. In the example considered, the isolated degenera- 
cies become a two dimensional degenerate subspace in 
which non-abelian holonomies are designed for adiabatic 
quantum computation. 



II. TOPOLOGY OF THE PARAMETERIZED 
EIGENVECTORS SPACE 

In this section, the topological features in parameter 
spaces are explained in simple words. Let a quantum 
system be dependent of N parameters x** defining a pa- 
rameter space P. Then, the Hamiltonian governing the 
dynamics is written as ?1I(R), where R = {x^ ,x'^ , . . . , x^) 
is a vector in P. The parameters are classical and can 
be tuned by an observer. Modifying the parameters 
amounts to trace a path F in the parameter space, pa- 
rameterized by time. To each point R in P is assigned 
the set {|/3(R))} of eigenstates of ?i(R), with the di- 
mensionality of the Hilbert space £. More precisely, for 
a single non-degenerate level a, a fiber J^^-' attach the 
eigenvector |a(R)) to the point R as: 

JT^' = such that 7i:(R)|a) = Ea(R)\a)^ . (1) 

This fiber is defined everywhere in P, except where the 
band a is degenerate. The set of all fibers attached to P 
defines the vector bundle T'-"^ over the parameter space: 
for a non-degenerate band a, it is a complex line bundle. 
A connection is a differentiable rule for a shift from the 
fiber to J^'^j^^-^ when R moves to R-l-dR in parame- 
ter space. When considering the adiabatic evolution of a 
quantum state |a(R)), parallel transport connections are 
involved. They arc such that 3m (c5;(R)|(iQ;(R)) = ev- 
erywhere along the path F covered. This requires that the 
path F never crosses a point where the band a is degen- 
erate. For a smooth choice of normalized states |a(R)), 
the parallel transport condition on a state \a) = e^'^"\a) 
is equivalent to a time evolution of the phase -ja = i{a\a) 
which can be integrated along the path F starting from 
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Rj as 



7a (t) 



R(t) 



(a(R)|Va(R)) dK. 



(2) 



This phase is purely geometric, i.e. independent of a 
reparameterization of coordinates on the path F. On the 
other hand, it is not invariant under a local gauge change: 
\a) — !■ e*^^^-* \a), making the phase 7^ non integrable and 
multi-valued. For this reason, the gauge field and its 
Berry's connection are specified as 



Af">(R) = (a(R)iVQ!(R)), 
= {a\da) = Af"'(R) 



•rfR. 



(3) 
(4) 



A key feature of the phase 7a is that it becomes a gauge 
invariant quantity when the paths F*^ are closed. In this 

case, Berry's phase^^ 



mod [2-jt] 



(5) 



is a physically observable quantity and cannot be re- 
moved by any local gauge change. It is sensitive to the 
topology of the fibre bundle: it is the holonomy^^ of the 
line bundle !F'-°'^ over the path F'^. An important special 
case arises when the Hamiltonian is real for a set of paths 
in a subspace of P: a continuous choice of real eigenstates 
|a(R)) may be chosen over this path, which defines a par- 
allel transport since 3to {a\da) = 0, leading to values of 
or TT for Berry's phase, such that |a(Rf)) = ±|Q;(Ri)). 

Berry's phase first appeared as the geometric contri- 
bution to the phase acquired in the adiabatic cyclic evo- 
lution of a non-degenerate state in the parameter space. 
In the next section, the state \ipa{t)) of a system initially 
prepared in the non-degenerate state |a(Ri)) evolving 
adiabatically along the path F is shown to be approxi- 
matively 



-»'/a(t)g«7c.(t) 



a(R(t))) , 



(6) 



where ria{t) = j- Ea{t')dt' is the usual dynamical 
phase and |a(R)), the instantaneous cigcnstatc. 

The gauge field A<°' defined in Eq. 4 is analogous to 
the vector potential of electromagnetism: in three di- 
mensions, the gauge insensitive magnetic field B*°' = 
V X A'°' is physically relevant. To characterize the prop- 
erties of the fiber bundle in a gauge independent man- 
ner, it is useful to define Berry's curvature as the dif- 
ferential form = B^^'J dx^^dx'^ = i{da\ A \da), where 
Bj^J = dfj_A\^^ — di^A'j^' arc elements of the antisymmet- 
ric curvature tensor . Using these definitions, Stokes 
theorem can be used to write the Berry's phase of the 
band a over the closed path F'' as a surface integral 



7a(r^) 



A'' 



(7) 



most pertinent for TQCP. The orientation of P is de- 
fined by a local choice of basis, for instance the natural 
basis (ui, U2, U3) with respect to the coordinates x^. Us- 
ing ordinary vector calculus, the antisymmetric curvature 
tensor reduces to a magnetic field B'"', the curl of 
A^"'. B<°' can be computed directly from the Hamilto- 
nian gradient as^^ 



\vn\a) 



(8) 



where {|/3)} is the set of eigenstates of H, dependant on 
R. As we shall see below, nonzero Berry's phases oc- 
cur from a non-trivial topology of the eigenvectors fiber 
bundle. This occurs at level degeneracies in the param- 
eter space where the magnetic field B<°' is singular. At 
these points V • B'"' ^ and A'"' cannot be defined. 
Without such points, the topology is trivial and the par- 
allel transport leaves states invariant over a closed loop: 
|a(Rf)) = |a(Ri)). 

The von Neumann- Wigner thcorcm-'^^ asserts that in a 
3D parameter space, accidental degeneracies may occur 
between two neighboring levels (say |±)) only at isolated 
points Ri*; these degeneracies which arc singularities of 
the fields B<*> have been named normal singular points 
by Simon^^. Since the gauge fields A<='=' do not exist at 
these points, they become local quantities. Interesting 
physics appear when such points live in P. It can be 
visualized most easily by projecting the Hamiltonian on 
the two- level manifold |±) which becomes degenerate at 
the singular point R* = {xl,xl,xl). Using the projector 
P = + ^ gradient expansion of the two- 

level projection H± of the Hamiltonian can be made in 
the vicinity of the singularity R* 



p[7i:(r)-?t;(r*)]P: 



H±(R) -H±(R*) 
Vn±{Il*)-SR + 0{dR^), (9) 



where SR = R — R* = {Sx^,Sx^, Sx-^). With a suitable 
choice for the origin of energies (£^±(R*) = 0), ?i±(R*) 
is zero. With this choice, this expansion can be expressed 
on the basis of Pauli matrices as 



n±{R) 



(10) 



The cf. are the elements of a 3 x 3 real matrix C which 
has a nonzero determinant for linear level crossing at R* . 
This becomes more familiar by defining b = Ci5R = 
(6^, b^jb^) as the effective magnetic field for an equivalent 
spin-i spin system, 



H±{R) = -o--b(R) 



b' 6^ - iby 
6^ + iby -b" 



(11) 



where S is an oriented surface with F'^ as a border. Sys- 
tems depending on a set of 3 parameters x'^,x'^,x^ are 



which magnitude increases linearly with the deviation 
from the degeneracy point b* = 0. The energy levels 
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E± = ±^ oi the two bands intersect conically at the 
degeneracy point (also called conical point or diabolical 
point). The matrix C maps a local neighborhood of R* 
of the parameter space onto a spatially isotropic spin- 
i hamiltonian in the magnetic field b. As long as the 
mapping amounts to a local deformation of the parame- 
ter space, and no additional degeneracies appear in the 
vicinity of R*, the topology of the fiber bundle stays 
unchanged. If C changes the orientation of space (the 
det(C') is negative), the sign of the topological charge 
is fiipped by the mapping. The one-to-one mapping C 
allows to use the Euler angles of b rather than the coor- 
dinate (5R to specify the eigenstates in the vicinity of the 
singularity as 



+ (b)) = 



cos : 



e"^ sm ■ 



■'o M-(b)) = 



-sm^ 
e"^ cos \ 



■ (12) 



For both levels, one can assign Berry's gauge potentials 



(13) 



These are the azimuthal gauge fields of "Dirac 
monopoles" of strength +i or — i placed at the origin. 
They are singular on their Dirac string {9 = t: for |-|-) 
and 9 = for |— )). These monopoles produce a radial 
magnetic field of opposite directions 



B(±' = Vb X A' 



(±) 



2|b|3 



(14) 



The strengths ± ^ are more easily identified by taking the 

divergence 



(15) 



which integrated over any volume including the origin 
gives ±i X 47r. Since the one to one mapping C between 
parameter space and spin-space conserves the fiux, the 
topological charge in P is preserved up to a sign (when 
C changes the surfaces orientations). 

The degeneracies (Ri*) in parameter space appear as 
singularities of the fields A'='=' . Since any surface S en- 
closing R* intersects the Dirac string, A<*' is not defined 
everywhere on S. It is possible to make A<*' single val- 
ued only by making a hole in S through which the Dirac 
string can be threaded: in this case, the surface can be 
continuously contracted to a point without crossing the 
singularity R* . This is the reason why there is no single 
analytic expression of A<*' over a surface which encloses 
completely the degeneracy. An alternative procedure for 
defining A<*> was made by Wu and Yang^'''. The space 
is divided in north (iV) and south (S) halves (see Fig. 1), 
with a diff'erent gauge choice A^='=' in each part, which 
are related by an appropriate clutching function on 
the equator where the eigenstates are connected using 




FIG. 1: A monopolo placed at the degeneracy points gener- 
ates azimuthal vector potentials A**' for each quantum band 
|±) intersecting at the monopolc. Because of the singularity, 
a single valued expression for A'^'^' exists separately in the up- 
per and lower half hemispheres. They can be connected on the 
equator by a clutching function (A'+' is pictured here). The 
corresponding magnetic fields B'*' are radial and decrease as 
b^^, where b is the distance to the singularity. 



|±)iv = e'^' '\±)s, with A^' = A^« + V/(±'. For 
the isotropic spin-i model, the diff'erent determination 
of A<*> are: 



A(±' = ± i 



cos 9 



2|b|sin6l 

cos 6* -h 1 
21 Y| sine 



— ^j\f 



for 9€ 



A(±) 



for 6'e (|;7r 



(16) 



with f-^^ = =F <;i as clutching function. 

On a closed path F^, Berry's phase for the two- levels ± 
is sensitive to the presence of a degeneracy at the origin 
since 



7±(n = £A<±'.d£ = T^f^(n. 



(17) 



where r2(r'^) is the solid angle seen from the origin. When 
F'' is contained in a plane intersecting the origin then 
7±(F'^) is just equal to tt times the winding number of F*^ 
around the origin. This discussion makes it clear that it 
is a consequence of the nontrivial topology of the bundles 
JF**' around the origin. A geometrical illustration is pos- 
sible when the path C shown in Fig. 2 lies into a plane 
where the Hamiltonian is real: the spin-eigenstates |±), 
which can be taken as real, depends on a single angle vari- 
able (say 9) which defines a line. As one moves along F° 
this line, which represents the eigenvector bundle, covers 
a one-twist Moebius strip as illustrated in Fig. 2. In this 
parallel transport, initial and final states are seen to be 
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l^(I*<)> 




paraUei^ transport 

FIG. 2; Illustration of a circular path C around a degeneracy, 
lying in a plane where the Hamiltonian is real. The real eigen- 
states depend on a single angle variable. This angle specifies 
the direction of a line. Under parallel transport along loop, 
this line (the eigenvector) generates a Moebius strip which 
is the visual representation of the eigenvector bundle. The 
eigenstate changes sign when coming back to the same point. 
This sign is the holonomy of the eigenstate fiber. 



opposite I ± (bf)) = — I ± (bi)): Berry's phase equals to 
TT. It is also the well-known property of the group SU{2), 
where rotations are Aw periodic. In general, if the path 
encircles m degeneracies and has a winding number 
around the degeneracy at the points R*, Berry's phase 
is mod [2tt] if X^I^i is even and n otherwise. 

Berry's phase factors for loops in planes where the 
Hamiltonian is real are topological invariants ±1, which 
characterize the sum of winding numbers around degen- 
eracies. Another topological invariant is obtained after 
integrating the field 6'=*=' of Eq. 14 over a small sphere 

around the origin. When normalized to 27r, these flux 
gives the first Chern numbers (or Chern indices) of the 
sphere with respect to the two bands |±): 

c[^\S^) = ^S B<±)-nd5 = Tl, (18) 



27r 



where h is the unit vector normal to the surface of the 
sphere. The Chern index can also be computed from the 
potentials A<*' using 



= T 



2n 



T^(27r-0)=Tl. (19) 



The mapping C between the parameter space P and the 
isotropic space does not change Chern indices if space 
orientation is preserved (det(C') > 0). 

Topological indices do not depend on the projection on 
a two-level system, which is valid only in a small neigh- 
borhood of R* . For any band a and any closed surface 



S in P, Gauss theorem assures that the integral of B'°' 
over the entire surface is identical to the sum of the inte- 
grals over a small sphere about each degeneracy. In other 
words, 



(<5)= E 9> 

diev{S) 



(a.) 



(20) 



where 17^'°' represents the topological charge of each de- 
generacy Rj* inside the volume V(S). The Chern index 
does not depend on the geometry of the closed surface 
and is a topological invariant which depends only on the 
degeneracies it contains. 



III. TOPOLOGICAL QUANTIZATION BY 
CONTROLLED PATHS 



Suppose that a quantum system depends on three tun- 
able parameters x^, and .t"^ which specify the space 
P. One can always construct two angles ■(9^(R), t9^(R) € 
[0, 2tt) in P which parameterize a two-dimensional torus 
and 1? ~ (i?^,i?^) is a vector on T^. Q is a physical 
quantity which can be expressed as the partial deriva- 
tive of the Hamiltonian ?i(i?) with respect to one of the 
angles (say ^?^) 



(21) 



TQCP can only be used for such physical observable, 
which is followed adiabatically on a path Fq lying on the 
torus. Physically, it is the ground state expectation value 
of Q which is of interest. On the path Fq, there will be 
one or more avoided level crossings with other levels, and 
Zener tunnelling in their vicinity sets the rates of varia- 
tion for the parameters required for adiabaticity. Using 
the spin representation (Eq. 11) close to a level crossing, 
where bz{t) is the tuning parameter, the condition for 
adiabaticity^^ may be written as 



(22) 



with a Landau-Zener transition probability Pl-z = 

exp ( — ^ |VJ ) . When this condition is verified for all 

V ^ hb^(t)J 

avoided level crossings on the path Fq, the adiabatic theo- 
rem^ may be applied to the non-degenerate state a which 
time evolution is approximatively 



IV'a(t)) 



aim)) , 



(23) 



where |Q:(i9(f))) is an instantaneous eigenstate (H|q:) = 
Ea\a)), the phase ria{t) is the usual dynamical phase 
factor 



Va{t) = I I E^{t')dt', 



(24) 



and ja(t), Berry's geometrical phase, was introduced in 
last section. In realistic systems, relaxation processes 
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restrict the use of TQCP to the ground state, and inelas- 
tic transitions to the first excited state will be shown to 
dominate quantization errors. 

Let us introduce the family of paths {T\} on dif- 
fering from Tq by a shift of t?^ by a constant angle 
A e [0, 27r). In the next sections, the helical family 

Tx-.tG [0; T] -dxit) ^ (27r i^i t, 2TiV2t + \) (25) 

where the angles and rotate at frequencies Vi and 
V2, will be used in a practical implementation of TQCP. 

When the frequencies are commensurate, the paths {r^} 
are closed. One of them is represented pictorially in 
Fig. 4, having commensurate frequencies z/i = 5 f2. Each 

path T\ begins at the point -di = (0, A) and the whole 
family {V\\ covers entirely the torus as A is swept from 
to 27r. 

The quantity of interest is the value of "the transferred 
charge" Q{T\) accumulated over the path {r^} which is 
covered in a period T. Integrating Q over time gives: 

Q{t)= [\d2nm')))dt'. (26) 

Jo 

The integrand in Eq. 26 is split in two parts 



(27) 



Each term contributes to the transferred charge Q{t): the 
first one leads to a dynamical contribution Q'^^", while 
the second one specifies the geometrical pumped charge 
Qgeo ifigntify the dynamical contribution, we take the 
time derivative of the adiabatic evolution (Eq. 23) and 
apply Schroedinger equation ihlipa) = WjV'a) to express 



j(?)a-7a) 



ih\a) + {Ea - h'ya)\a) . (28) 



The expectation value of H and its phase derivative fol- 
low from Eq. 28 and Eq. 23, 



d2{H) = d2E^ + nd2{i{a\a) - ia) . 



(29) 



Since 7q = i{a\a), the last two terms on the right hand 
side disappear. When integrated over the period T this 
first contribution Q'*^"(rx) to the transferred charge is 
also the derivative of the dynamical phase with respect 
to the initial angle A 



g'^^"(r> 



= rd2E^{^^{t))dt=h^'^''^^^^ 

Jo 



dX 



(30) 



where the definition (Eq. 25) of the helical paths has been 
used to transform the partial derivative of the integrand 
into a total derivative of the dynamical phase with re- 
spect to the initial angle A. This quantity is just the dif- 
ference between the total accumulated dynamical phases 
on the neighboring paths F^v+d;^ and Tx normalized to 
the angle increment dX. 



We now turn to the geometrical contribution Q^^°, 
which comes from the second term in Eq. 27. Taking 
the 1?^ derivative of Eq. 23 yields 



\d2'ipa) 



\d2a) + id2{'^a - Va)\(^) ■ (31) 



Using Eqs. 29 and 31, the scalar product {tpa\'H\d2'>pa) 
gives several terms, but only one of them is not purely 
imaginary, namely —ih{ce\d2a). When integrated over 
time, the result does not depend on the dynamics, but 
only on the path geometry. Hence the geometric pumped 
charge Qf°{rx) is 



Q«-(rA) 



-2 3? / {Mt)\n^xmd2Mt))dt 
Jo 

h [ 23(a2Q!(t?(i))|rfa(i?(t))). 



(32) 



This charge can be expressed in term of a geometrical 
phase by rewriting 



2Qm {d2a\da) = i d{a\d2a) — d2{a\da) 



(33) 



The second term is recognized the derivative of the 
connexion Aa (defined in Eq. 4) which integral over a 
closed path is Berry's phase. When integrating over the 
path Tx, the first term only contributes at the endpoints 
i?i = ■i?a(0) and 'df = 'dx{T), giving 



Qg-(rA) 



Jtx 



ih {a{-df)\d2a{-df)) 
- {a{^i)\d2a{^i))] , (34) 



these last two contributions being essential to enforce 

the gauge invariance of Q^'^°{rx), a measurable quantity. 
When the path is closed, the endpoints contributions 
cancel, and Q^'^°{Tx) is the integral of the t?^ derivative 
of the vector potential 



Qs<^°{Tl) = -hd)d2A 



1(a) 



(35) 



For the helical path family {F^}, we showed in Eq. 30 
how the dynamical transferred charge Q'^^'^{Tx) could be 
expressed as the total derivative of the dynamical phase 
with respect to the initial angle A. The same argument 
can be used here mutatis mutandis to the geometrical 
transferred charge 



-^Tx 



A'-"^ = -h 



dX 



(36) 



This formula presents the advantage to be easily gener- 
alized to open paths thanks to the endpoint contributions 
in Eq. 34. In a first step, the integral of the angle deriva- 
tive 



A 



(37) 



FIG. 3: The virtual path is the sum of the four path- 
segments T\, Tf, and ^ lying on the torus The 
geometrical contribution Q^"" to the charge Q over the seg- 
ment Fa is proportional to the circulation of A^"' along the 
"virtual" closed path S^. 



FIG. 4: Representation of a closed helical path F^^ lying on 
the torus (here with a frequency ratio 1^2/1^1 = 5). The 
geometrical contribution Q^"" of the charge Q over a period 
T2 = ^ is proportional to the flux of B*°' through the in- 

finitesimally thin strip S^^, . 



is rewritten as a difference between two paths shifted 
by the infinitesimal i5A, which becomes a sum when one 
of the segment is integrated in the opposite direction 
(r^^^_^). These two paths can be connected by infinites- 
imal vertical segments and F/ at their endpoints i?i 
and 1?/ as shown in Fig. 3. The endpoints contributions 
in Eq. 34 can be rewritten as line integral of the vector 
potential over these end-segments as 

(38) 

-ih{a{^i)\d2a{^i)) = 

When combining the four path-segments together, a 
closed path is constructed from the path V\ which 
is one of its line-segments as drawn in Fig. 3. On this 
closed path S;^, the integral of the vector potential be- 
comes precisely Berry's phase. By constructing the four 
segments virtual path S^, one of which is the physical 
path F;^ of interest, the geometrical transferred charge 
on Fa can be written as 

Qg-(FA) =ft^^2g>^. (39) 

In contrast with Eq. 36 which gives only a global de- 
scription of the geometric charge on a closed path, this 
expression for the pumped charge can be used on any 
arbitrary paths. They are relevant if noise or error in the 
control of parameters exist. 

What is the benefit of this formulation in term of 
Berry's phase ? One is practical: Berry's phase can be 
computed efficiently. The gauge dependence of the vector 



potential A<°' introduces a difficulty which can be cir- 
cumvented in two ways. Berry's phase can be computed 
as the flux of the magnetic induction B*°' using Eq. 7. 
The two dimensional integration can however be tedious 
to compute, particularly when the surface is warped. Al- 
ternatively, King-Smith, Vandcrbilt and Resta-*^^'^" for- 
mulated Berry's phase in term of a gauge invariant ex- 
pression by discretizing the one- dimensional path Fa. Let 
'&j be N points on Ea splitting it in A'' small segments. 
The line integral of the vector potential can then be ex- 
pressed as the invariant^^ 

JV-l 

7a(SA) « - arg n (a(i?,-)|a(^,+i)) • (40) 

In this way, any local gauge change cancels out between 
bras and kets which come each in pairs. It is easy to 
implement over complex paths and very accurate. When 
shrinking the path Fa to an infinitesimal segment be- 
tween 'd\{t) and 'd\{t + dt), one also gets the instanta- 
neous geometrical pumped charge: 

<5g«°°(i) ^ argn(a(t?,)|a(^,+i)> , (41) 
3=0 

where the four points are the extremities of the infinites- 
imal paths (sec Fig. 3) . The geometrical pumped charge 
Q, which is a physical quantity, can be tracked and mea- 
sured anywhere along any real path. Equations 39 and 41 

are thus of great practical value since the local physical 
processes and experimental sources of errors in the path 
can be analyzed on the quantity of interest Q^°°. 

For particular sets of paths (the "controlled paths"), 
the "charge" Q = (SaTi) can be quantized through its 
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relation to the Chern index of a closed surface for the 
ground state eigenvector bundle. Any closed surface 
containing one or more singularities can in principle be 
used. For simplicity, the entire torus will be used 
here. There are several ways one can generate this two- 
dimensional surface using a one dimensional path. The 
helical family {Tx}, is one of the possible families of 
controlled paths which generates the surface T^. When 
the angular frequencies of and arc commensurate 
1/2 = pvi, the angle winds p times around in a pe- 
riod. When p is large, the helix covers densely the torus. 
Alternatively, the initial angle A can be swept from to 
27r to sweep the helix on the torus surface p-times. Using 
this averaging procedure, the dynamical contribution to 
the pumped charge averages out to zero. 



dX 



dX 



(42) 



/l(r?„(r^J-ry„(rg)) =0, 



since = Fq. When discussing the geometrical con- 
tribution (Q'^''"), it is simpler, to split the helix into p 
one-turn segments r]^ (which are open paths). As 
the initial angle A of the helix is swept from to 2tt, each 
one-turn segment F^ sweeps the torus surface just 
once {p times for the whole helix). For this reason, it is 
simplest to compute the average of the geometric charge 
over a 2tt A-period for this one turn segment F^ and 
multiply the result by p for the whole helix. The integral 
over A can be made by dividing the 27r period in A'' small 
slices indexed by j of width AA = 

The contribution to the pumped charge over the one 
turn segment rj^^*'"" for the j**^ slice of width AA defines 
the helix strip of surface Sl^ . represented in Fig. 4. Its 
boundary is nothing but the virtual path Sj '"'''^ associ- 
ated to the one- turn segment Fj^^*"''" (see Fig. 3). Using 
Eq. 39, the pumped charge averaged over this interval is 



A/^gco/T^l turn\ „, /v-'l turn\ 



(43) 



/ A<°> • dR 

Z\A Jj2\ turn 



where Stokes theorem was used to express Berry's phase 
along the virtual path as the flux of B^"' through 

. When summing over all the j slices of height AA, 
these elementary surfaces add up to the entire surface of 
the torus. Hence, when averaged over A, the geometrical 
charge transferred becomes 



N 



N 



^2:^"(^^-) = ^FiB--»*5^ (44) 



This last term is precisely the Chern index c^°' of the 
surface for the a-eigenvector bundle. For the whole 
helix, each one-turn segment contributes equally and 



(45) 



The average of Qseo ^y^j. ^j-^g family {F^} is quantized by 
the winding number p, and the Chern index of the torus 
with respect to the band a. It is nonzero only if degen- 
eracies involving the band a are present inside the torus. 
As was pointed out by Goryo and Kohmoto^^, invariances 
of the Hamiltonian under mirror symmetries (i?i,i?2) 
{{-^1,^2) or i^i,-^2)} (i.e. n{^i,^2) = ni-^u^2) 
or 'H{'di, —'02)) and "time-reversal" symmetry 1? — 1? 
(i.e. 'H{'&) = H* {—'&)) are incompatible with a nonzero 
Chern index. This is because B<°' is an axial vector: the 
mirror symmetry (z?i,z?2) (^^1,-^^2) leaves the torus 
invariant, but B'"' changes sign with respect to the local 
natural basis (e^i , e^2 , n) (n is the vector normal to the 
surface) . Hence the mirror symmetry switches the sign of 
B<°' (i?-'^, i?^) • n and the integral of B<°' • n over the torus 
vanishes. This property can be used locally to detect the 
presence of singularities in the eigenvector bundle. For 
example, in the spin representation (Eq. 11) close to a sin- 
gularity, under the mirror symmetry {bz —bz) the |-|-) 
and |— ) are mapped into each other and each eigenvector 
bundle is not preserved separately. The same behavior 
occurs under time-reversal. 

Goryo and Kohmoto^^ generalized the relation be- 
tween the expectation value of a derivative of the Hamil- 
tonian and the Chcrn indices on D-dimensional torii, with 
application to a number of problems (IQHE in 2D and 
3D dimensions, ACJE, etc.). In these problems, the av- 
eraging over the whole torus can be made directly, but in 
our case the physical quantity Q is generated by paths: 
a path description cannot be avoided. Since the average 
(Q) over a family of commensurate paths {F^} is quan- 
tized, the value of Q(r^) for a given A fluctuates around 
the integer mean value. It is interesting to know how 
these fluctuations decrease with winding number. Since 
the torus is covered densely at large p, we expect a more 
acciiratc quantization as the winding number p get larger, 
irrespective of the value of A. A more accurate averaging 
of the dynamical charge improves the quantization. This 
will be easiest if t?^-dependence of the energy Ea{i3'^) is 
weak since Q^^^" ~ d2Ea ■ TQCP is an asymptotic quan- 
tization, which works best for the ground state which is 
most robust against incoherent processes. For a two- 
dimensional torus, a number of paths can be chosen, the 
only requirement for TQCP being the 1?^ periodicity. 

Next section, devoted to the Cooper Pair Pump is a 
physical example where TQCP can be implemented con- 
cretely. 
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FIG. 5: The circuit is made of a Cooper Pairs Pump closed 

on a small inductance L, threaded by a magnetic flux $. The 
charge on the islands 1 and 2 can be tuned by two gate volt- 
ages Vgi and Vg2. The inductance is used to bias the phase 
across the CPP, </? = 27r$/^'o. 



IV. TOPOLOGICAL PROPERTIES OF THE 
COOPER PAIR PUMP 

One of the simplest implementation for a Cooper pair 
pump (CPP) using a superconducting circuit is repre- 
sented on Fig. 5. Phase biasing is achieved by closing 
the CPP on a small inductance L, threaded by a mag- 
netic flux <I>. Its magnetic contribution to the energy is 
- ip)"^ where (j) is the phase difference across L. For 



1 

2L 



small L, it has a deep minimum at (p — 2tt-^: this in- 
ductance and the CPP series capacitance Cg form an har- 
monic oscillator which frequency {2tt^/LCs)~^ exceeds all 
other energies, effectively blocking the quantum variable 
(j) at the value (the center of the ground state wavefunc- 
tion). (p is then a parameter tunable by the magnetic flux 

The three Josephson junctions, with small capaci- 
tances define two superconducting islands with suSicient 
large electrostatic energies to limit charge fluctuations 
through the junctions. Let be the excess number of 
Cooper pairs (with respect to charge neutrality) on is- 
land i. The electrostatic energies of each island can be 
tuned independently using a gate voltage Vg^ through 
the gate capacitances Cg. . The induced charge polariza- 
tion on the island i is ng. = Cg.VgJ {2e) in units of 2e. 
For convenience, we use the total charge n., = ni + n2 
on the double-island and the charge asymmetry Ud = 
n\ — 122 between them as the natural basis of charge 
states {[rig, na)}. Taking the two external junctions with 
the same Josephson energy Ej and capacitance Cj, and 
{EjQ, Co) for the central junction Josephson energy and 
capacitance, the charging energy of the CPP reads 



He = Ec 



{fls - 



(46) 



where Eq 



{2ef 



Ugd-- 



'91 



n, 



iCj 

92 

1 ^ 



ga- 



''91 



''92 ) 



is the Coulomb energy, n, 
and Ko = 2C^+Cj ^ capacitance ra- 



tio (of order -i). In addition to the phase bias p, the 
induced charge ngs,ngd are tunable parameters of the 



Hamiltonian: the parameter space P is here three dimen- 
sional, and a point R in P is specified by its coordinates 
{ngs,ngd.,p). One easily checks that the charge state 
\na, Ud) is the ground state which minimizes the parabo- 
las in TLc (Eq. 46) inside the hexagonal area h(^ns,na) cen- 
tered at the point {ngs= ns,ngd= rid), in the Ugg — Ugd 
plane (Fig. 6). On the line boundaries between hexagons, 
two electrostatic states have the same energies while the 
vertices are points of triple degeneracies. This hexag- 
onal lattice of triple degeneracies has two points in its 
unit cell, chosen here as {Ti} — {rigs = ^-^,ngd = O) 
and {T2} = {rigs = ^~2° , rigd = l)- The Josephson tun- 
nelling, which "translates" Cooper pairs across the junc- 
tions, can be expressed in term of the variables conjugate 
to the total charge and charge asymmetry rid, ©s and 
Qd which are the generators of charge translations, 

Hj = -2Ej cosQsCos{Qd- Kj p) - Ejq cos{2Qd + kq (p) , 

(47) 

where kj = 2Co+Cj ~ other capacitance ratio 

(also of order i). Since Hj delocalizes Cooper pairs, the 
charge states are no longer eigenstates of the full Hamil- 
tonian H = He + Hj and the degeneracies along the 




FIG. 6: Stability diagram for the charging Hamiltonian 
(Eq. 46) in the rigs — rigd plane. The charging energy is 
minimized inside the hexagonal areas shown. The 

boundaries between hexagons are line of degeneracies between 
charge states, while the vertices are points of triple degenera- 
cies. The two vertices Ti and T2 form the unit cell for this 
hexagonal lattice. The coordinates of Ti are ( , 0). Also 
shown are the topological charges of the lattice of degeneracies 
in the plane p = tt (see text). 
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boundaries of the honeycomb lattice are Hfted. Neverthe- 
less, if Ej « EjQ < Ec, accidental isolated degeneracies 
persist in P in the vicinity of the points Ti = {Ti, (p — tt} 
and T2 = {T2,'!t} and all their equivalents under lattice 
translations in P. In the special case where Ej = Ejq 
(homogenous array), the degeneracies are placed at Ti 
and T2, and shift along the rigs-axis for asymmetric ar- 
rays. As an illustration, the energy manifolds for the 
two lowest levels |±), represented in Fig 7, show the two 
conical intersections in the points Ti and T2 in P for 
Ej ~ Ejq. Symmetries of the total Hamiltonian are 
most explicit after the unitary transformation generated 



by Uiip) 



' , which leaves the charging Hamilto- 



nian and shifts Tlj into 

nj^u{ip)HjU\cp) 

= -2Ej cos Qs cos Qd - Ejq cos(2ed -1- ip). (48) 

In this representation, the phase bias appears across 
the central junction instead of being distributed across 
the three junctions according to the electrostatic volt- 
age drop. The mirror symmetry n 



keeps the 

physics unchanged while {hs,6s) {—ns,—9s). When 
the phase f is equal to or tt, the Hamiltonian is 
real and also invariant under a second mirror symmetry 
rigd -rigd, while {nd,Qd) {-fid.-Qd)- The inte- 
ger translations on the honeycomb lattice (ngs,ngd) — > 
{ugs + UsQ, Hgd + ndo) induced by the translation opera- 
tor e"*'^"'"®""''"'*"®''^ leads to different numbers of Cooper 
pairs on the island {ug — Uso, rid — rido)- These states are 
equivalent but physically distinguishable. Finally, the 
Hamiltonian is 27r-pcriodic in ip and phases differing by 
multiple of 2tt leads to identical physical states: ip plays 
the same role here as the i?^ variable in the preceding 
section. 

In the same fashion as in section II, let us construct 
explicitly the two-levels approximation of Ti. in the vicin- 
ity of the degeneracy point Ti. Since circuit asymme- 
tries do not affect the topology of the eigenvector bun- 
dle, it is simpler to take symmetric junctions Ej — Ejq 
and assume that the ratio P = between Joseph- 

son and charging energies remains small. Writing the 
small deviations from the triple point Ti as u = Ugs — |, 
S = Ugd, Ip = p — TT, the projection of the Hamiltonian 
ii on the basis of charge states {|0, 0)', |1, 1)', |1, — 1)'}, 
{\ns,nd)' ~ U(p)\ns,nd)), is represented by the matrix 



n^E, 




-fi -(3 
-2^^ (3(1 -i^) 
(3 (3{l + iTp) -2^ 



,{a-5)_ 

3 



to first order in the deviation (5R = {a,6,ip). 
{a = S = -i/j = 0), the two lowest eigenstates 



(49) 



At Ti 



1 



1,1)'- 



V3 + 1 
2 

Vs-1 



1,-1 



1,-1 



111 
111 




FIG. 7: (a): The two lowest energy manifolds computed in 
the rectangular area of the rigs — rigd plane shown in Fig. 6 
for the constant phase p — tt. In the vicinity of the isolated 
degeneracies Ti and T2, the energy sheets form a conical 
intersection, referred to as "diabolical points" (same shape as 
a diabolo). At these points there are topological singularities 
in the bundle JT'-' (Eq. 1). 



are degenerate with energy — (ground states) , and the 
first excited state 



k(T,)) 



^(10.0)' 



11,1) 



11,-1) 



(50) 



has Ej for eigenvalue. As discussed in Sec. II, an 
isotropic spin representation of the Hamiltonian in the 
|±) subspace require a deformation of the parameter 
space P represented by the matrix. 



C 













~^E 








(51) 



which amounts here to a symmetry (the flip of the S 
and ■(/' axes changes the space orientation) and a linear 
deformation. This transformation specifies the effective 
magnetic field (b = CSR), b'^ = ^Eca,by = -^Eji/j 

and b^ = |^ S, such that the projection of the Hamilto- 
nian on the |±) degenerate subspace reduces to a spin-i 
Hamiltonian 



n 



±(Ti 



&^ 6^ - i by 

b"" + iby 



1 



cr-b(R). (52) 



The two lowest levels have a conical intersection at the 
degeneracy point Ti 



£;±(b) = ±y 



(53) 
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FIG. 8: Comparison of position of the degeneracy point Ri* 
between the analytic model described in the text and the 
exact numerical result. 



Following the discussion in the Sec. II, the topological 
charge in the spin representation and in the original pa- 
rameter space are identical up to the sign of the determi- 
nant of C which is positive. Hence the topological charge 
of the ground state is (/'"'(Ti) = +1 and —1 for the first 
excited band. Using the same arguments, the topological 
charge at the degeneracy T2 is (j''°'(T2) = — 1. Similarly 
all degeneracies obtained by lattice translation from Ti 
(resp. T2) have a topological charge of +1 (resp. —1) for 
the ground state. As mentioned in Sec. II, the transfor- 
mation properties of the eigenstates bundle under mirror 
and time-reversal symmetries allow to detect the pres- 
ence of a degeneracy locally (in the spin-i representa- 
tion). Here, the |±) states map into each other, and 
each eigenvector bundle is not preserved separately by 
these transformations. When Ejq deviates from Ej, the 
degeneracies move continuously away from Ti and T2. 
Using the same 3x3 matrix representation for H, the 
degeneracies slide along the ngs axis in the vicinity of Ti 
and T2 as 



R* 
1,2 



(Ti)T 



1 



'iEc 



E 



JO 



El 



,0,7r 



(54) 



for small deviations \Eja — Ej\. In this shift, the degen- 
eracies keep their topological charge, (/'"'(R*) = g<'"(Ti). 
In Fig. 8, the analytic and exact positions of the degener- 
acy points are compared: the agreement deviates rapidly 
as one moves away from Ti . 

The presence of degeneracies in the lowest band of the 
CPP allows to quantize the pumped current opening ac- 
curate application for metrology, the topic of next sec- 
tion. 



uses normal island separated by tunnel junction were bi- 
ased by a small dc voltage Vb- If the two gates volt- 
ages are driven in quadrature, the systems undergoes a 
circular cycle centered around point Ti. The electro- 
static ground state changes cyclically (|0, 0) |1, — 1) ^ 
|1,1) |0, 0)) as one crosses one of the three degen- 
eracy lines intersecting at Ti. After one cycle, a single 
charge is transferred through the electron pump. If the 
cycles are sufficiently slow, the charge relaxation (e.g. 
|0,0) — » 1 1,-1)) is inelastic but has sufficient time to 
complete. Since the process is stochastic, errors occur 
and limit the accuracies of normal electron pumps. Also, 
the timcscale for charge relaxation are typically of order 
T = {RtC)~^, where Rt is the tunnelling resistance. For 
realistic circuits, r rarely exceeds 10~^ sec, and pumped 
currents do not exceed a few pico- Ampere. 

This is on of the motivation for studying Cooper-pair 
pumps^^"^^ (CPP) to circumvent the stochasticity of nor- 
mal electron devices. Here, we show that the charge 
transferred can be quantized topologically by using con- 
trolled paths in parameter space (TQCP). The CPP's cir- 
cuit delivers a current I which is equal to the charge 
transferred Q per unit of time : I = Q (see Fig. 5). Let's 
return to Eqs. 46 and 47 and consider that (p is still a 
quantum degree of freedom conjugated to a charge oper- 
ator q, i.e. [q, 0\ = i. The time evolution of the mean 

value of q is = —j[{[q,'H]), and is equal to ^ {d^H). 
Since the small inductance blocks the quantum fluctu- 
ation in (f, it can be taken as classic and the pumped 

current is 



/ = Q = f (d^n) . 



(55) 



Since if and 1?^ have the same 2tv periodicity j^I has the 
exact expression (Eq. 21) as required for the TQCP pro- 
cedure discussed in Sec. III. Consider now the cylinder 
<S in parameter space represented in Fig. 9, which axis 
lies in the direction. Its section in the Ugs-Ugd plane 
has a radius p of order ~ | and its liciglit on the 
<y9-axis is 27r. Since the end-faces p = and p — 2Tr are 
physically equivalent, this cylinder <S is a closed surface 
and has the topology of a torus . A point on S is spec- 
ified by two angles "ff^, the angle in the Ugg-Ugd plane, 
and the phase p = {F. The cylinder's radius p is cho- 
sen so as to include only one degeneracy Ti. It has the 
same Chern index as any other surface which includes 
Ti : cf{S) = g'^HTi) = +1 with respect to the ground 
level. By deformation, the cylinder offers the advantage 
to contain the same helical paths (Eq.25) as the one on 
the torus used in Sec. HI. TQCP can therefore be used 
exactly in the same fashion for the topological charge 
which is quantized as 



V. QUANTIZATION IN THE COOPER PAIR 
PUMP 

The idea of using single electron pumps^^ as ciirrent 
standard dates from the early 90's. The original circuits 



(56) 



when averaged over the initial phase A of the helix (de- 
fined in Eq. 25) making p turns around Ti. 
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FIG. 9: A helical closed path To in P on the surface of a 
cylinder enclosing the degeneracy Ti. Since (p is cyclic, the 
surface 5 is a two-dimensional torus and To is closed for in- 
teger frequency ratio p = ve/v^ (here 2). 



Since Eq. 41 gives the incremental charge transferred, 
this quantity can be monitored continuously as one moves 
along the helical open path using Eq. 39. The Josephson 
couplings Ej and Ejo ~ Ej lift the degeneracies be- 
tween the charged states along the three boundary Unes 
of the hexagons intersecting at Ti, which become saddle 
lines (see Fig. 6 and Fig. 10-b for the saddle lines of the 
ground state). The splitting between the ground state 
|— ) and the first excited state |+) is smallest along the 
saddle lines and are of the order of Ej. As one moves on 
the helical path, each saddle line is crossed once per turn 
around Ti . Since the accuracy of the CPP hinges on the 
ability to move adiabatically in the ground state mani- 
fold, Landau Zener- transition |— ) \+) when crossing 
a saddle line are a concern. The transition probability 

Pl-z — e ^ ^ ' ^c- depends on the ratios ^ and 
which cannot be too small. On the other hand, when 
the ratio ^ is too large, the dynamical contribution 
to the charge transferred are more difficult to average 
out and the accuracy of the device deteriorates. This is 
the tradeoff when optimizing the CPP: a large Ej re- 
duces Landau-Zencr tunnelling and allows for a higher 
frequency of operation but the Josephson current can be 
most easily driven to zero at small Ej. To avoid single 
electron effects, the charging energy Ec has to be smaller 
than the superconducting gap (0.2 meV for Ahuninum). 
This sets the overall energy scale and most of the param- 
eters: typically values for Eq « 0.1 meV, Ej « 0.05 i^c 
and U0 w 100 MHz offer a good optimization of the CPP. 
The parameters of the helical path are the radius p and 
the number of turns p in a period = |^ . For the op- 
timal radius p = « |, the path intersects the line 



between Ti and T2 in the middle. A smaller radius is 
equivalent to reducing Eq- 

Now that the parameters are known, the charge trans- 
ferred can be followed as one moves along the helical 
paths defined in Eq. 25. p is here the number of turns 
aroimd Ti in a period T^. Integrating the instantaneous 
transferred charge (Eq. 41) for the ground state |— ) yields 
the time dependence of Q^^° 

Q^^°{t)= f 5Q^''°{t'), (57) 

which is plotted as a function of time in Fig. 10 using 
the parameters of a real device (^ ~ 0.05) and Eq. 41. 

For small ^ ratios, the charge is transferred in three 
distinct steps, corresponding to the transfer of a Cooper 
pair through each junction which occurs when crossing 
the three saddle lines. Two steps are of height ~ kj 
(external junctions) while one is of height ~ kq, yielding 
a total transferred charge ~ 2kj + kq = 1 per turn (in 
units of 2e), as illustrated in Fig. 10 where this quantity 
is plotted for different number of turns {p). The steps 
rounding become more pronounced as Ej increases, and 
their size more sensitive to the initial phase value (p(0) = 
A. 

Clearly, the charge transferred value ~ 2e is due to the 
presence of the degeneracy with topological charge 
We now verify that quantization accuracy improves as p 
increases, a fundamental feature of TQCP: 

g(rA) 2ep. (58) 

p — ^00 

The geometric and dynamical charge transferred in the 
CPP ground state can be followed as a function of the 
initial phase A. For optimal values of the parameters 
{Ec « lOOpeV and Ej « 3peV), the errors computed 
are small. In the simulation, it is useful to amplify their 
effects by choosing the most unfavorable parameters. Us- 
ing a perturbative analysis, the dynamical contributions 

to the pumped charge are of order Q'^^'^ w in 

units of 2e^®. Similarly, the deviations of the geomet- 
rical pumped charge from its quantized value scale as 

Hence, larger values of Ej/Ec increase errors. In 
Fig. 11, the geometrical and the dynamical charge trans- 
ferred are plotted for a ratio Ej/Eq = 0.5, an order of 
magnitude larger than the optimal values. On this figure, 
the charges are computed using helical paths around Ti 
with different number of windings in a period T^. The 
geometrical charge oscillates as a function of A around 
the quantized value with an amplitude which decreases 
rapidly with the number of windings p. This rate de- 
pends mostly on Ej/Eq and on the distance between 
each windings relative to the distance to the degeneracy 
Ti. This is the main reason to keep the helix radius close 
to its optimal value { ^~2° )- 

For parameters closer to their optimal value, this de- 
crease can be expressed in term of ^, the root mean square 
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t/T 



FIG. 10: (a) Geometric transferred charge Q^°° as a function of t for the hehcal paths defined in Eq. (25) and p = 1, 3, 5. (b) 
Each time the path crosses one of the three saddle lines (here plotted for Lp = n), approximatively a third of one charge unit is 
geometrically transferred. The computation was done using a radius p = 0.3, A = 0, and the Josephson couplings and charging 



energies 


of a real device: Ej 


= 6 neV, Ejo = 


8.6 neV, Cj = 


P 


Cgco 


^dyn 


^tot 


1 


8.7 lO-" 


5.4 


5.4 


2 


2.8 10"^ 


4.2 10"^ 


5.0 10"^ 


3 


3.8 10"* 


1.3 10"^ 


1.3 10"^ 


4 


3.9 10"^ 


9.1 10"^ 


9.9 10"^ 


5 


4.6 10"^ 


5.5 10"® 


7.2 10"® 



TABLE I: Mean root squares deviations of the different contri- 
butions to the transferred charge as a function of p. Numeric 
values : p = 0.3, \ = 0, Ej = 60 p,eV , Ej/Ec « 0.5 . 



amplitude of the oscillations. This quantity is tabulated 
as in the table I using Ej/Ec = 0.05. Above a few wind- 
ing the quantization accuracy is very high. If low fre- 
quencies phase jitters in 6(p are present, the error in the 
pumped charge will be of the order of ^geo(p)^) which 



is below 10 ^ for p > 5. 



The average dynamical charge transferred over a pe- 
riod converge also toward zero when p is sufficiently large, 
provided the ratio of is not too large (say below than 
0.05). In absence of noise, the periodicity in guarantee 
that it averages out to zero. In the presence of a phase 
noise Sip the cancellation becomes approximate with an 
error of order ^geo(p) 



Sip 
27r- 



It is not possible to reduce the ratio Ej/Ec arbitrar- 
ily to improve the accuracy, because the gap at the sad- 
dle points (oc Ej) decreases, and the Laudau-Zener tun- 
nelling turn on transition to the first excited manifold. 
This introduces the largest source of errors because the 
Chern indices of the two lowest eigenvector bundles are 
opposite, an issue which is addressed in the concluding 
section. 
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A (bj A 

FIG. 11: (a) Plot of the geometric charge Q^"" as a function of of A for helices with different numbers (1-4) of windings, (b) 
Same plot for the dynamical charge Q'^^'^ as a function of A using 4 to 7 windings. Parameters: p = 0.3, Ej = Ejo = 60 //eV 
and Ec = 114 //eV. 



VI. ADIABATIC COMPUTING 

In the vicinity of a triple point Ti, the two lowest states 
I ± (R)) form a qubit. In this region of P, the next level 
|e(R)) Ues at an energy of order above the |±) dou- 
blet, significantly larger than the doublet splitting. One- 
qubit operations on this doublet can be implemented us- 
ing different schemes. 

The simplest one consists in applying microwave pulses 
to the CPP gates at the frequency splitting between the 
I ± (R)) states. The point R may be chosen at one of 
the magical points along the saddle lines of the hexagons 
where the system is to first order insensitive to gate 
charge and phase fluctuations^. But there are simpler cir- 
cuits where such operations have been demonstrated and 
the CPP is more interesting to implement quantum com- 
putation using adiabatic cycles. Duan and coworkcrs^^ 
have shown how geometric gates could be implemented 
on a degenerate two-levels system using resonant transi- 
tions between this |±) doublet and an excited level |e) 
provided it is also coupled to an auxiliary level \a). In 
the CPP, I a) could in principle be a higher lying charge 
state, but this solution is not as convenient as for atomic 
systems. We prefer to do without the \a) state. In this 



case, the dynamical contribution associated to the Rabi 
frequency cannot be eliminated altogether. As long as 
this phase shift can be tuned to a multiple of tt, a geo- 
metric gate can be implemented. 

The CPP is biased at degeneracy (Ti), and microwave 
voltages tuned at the |e) — |±) frequency splitting arc 
applied to both gates. When going to a frame rotating at 
the same frequency, the coupling to the |e) state becomes 
time independent 

n^^it) = m (^cos^|e)(-| +sin^e''^|e)(+| +h.c.^ , 

(59) 

where i7 is the main Rabi frequency, and the angle Q con- 
trols the relative strength of the couplings of the | -h) and 
|— ) to the |e) state, while x is their relative phase. These 
three quantities can be adjusted by tuning the amplitude 
and the phase of the microwave voltages applied to the 
CPP gates. Here, 2tt phase cycles in x are sufficient to 
generate the gate operations. In the rotating frame, one 
of the sate is stationary, while the other two oscillate at 
±0. Adjusting the period of operation T to eliminate 
the dynamical phase shift (e*^^ = ±1), a 27r-cycle in x 
generates the operation 



G(C,27r) = 



_ / cos(7rcos^) cos^ -|- i sin(7rcosC) 



cos(7r cos C) sin C, 



— cos(7r cosC) sin^ 
cos(7r cos C) cos C — * sin(7r cos Q) 



(60) 



which covers all the one-qubit gate operation. 

Two-qubit geometrical gates can also be considered by 
coupling two Cooper pair pumps together. One way this 

may be achieved is to close the two pumps on the same 
inductance L: the current pumped in both device add 



up in I/, shifting the phase ip by an amount of order Llj 
where Ij is the current through the j-th CPP. 

In practice, adiabatic computing faces a number of dif- 
ficulties. Even for one-qubit operations three or more 
states have to be degenerate. But any low frequency 
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flux or charge noise moves the area in parameter space 
where the degeneracy occurs and the shifts in parameters 
rotates the eigenvectors rapidly. While this has little im- 
pact for topological quantization, this deteriorates the 
performance of geometrical gates. The accuracy with 
which parameters are to be controlled is also consider- 
ably higher than for usual quantum gates. These are 
some of the reason why geometrical gates have not yet 
been demonstrated. 



VII. CONCLUSIONS 

Most limitations for adiabatic computing are irrele- 
vant to topological quantization which is robust against 
most perturbations. For a CPP, charge or phase noise 
is considered slow when most of its spectrum is below 
the frequency of operation of the pump (100 MHz is a 
realistic number). Such noise source add a random com- 
ponent to the controlled voltage or the magnetic flux and 
the path in P no longer generates a cylinder, but a more 
irregular surface. As far as the geometrical charge is con- 
cerned, this has almost no effect, as long as the resulting 
surface still encloses the topological point Ti. Similarly 
if the junction capacitances or Josephson couplings fluc- 
tuate in time, the position of the point Ti fluctuates in 
P, and this has little effect as long as this shift is small 
compared to the cylinder's radius. In presence of low fre- 
quency noise, dynamical contributions no longer cancel 
exactly. On the other hand, the errors are random and 
can be averaged using long integration times. 

If charge or phase noise has frequency components at 
the splitting between the two lowest states | ± (R)) or if 
Landau Zcner transition occur at one of the three saddle 
point crossings (see Sec. V for a discussion), the system 
will spend a fraction of the time in the first excited state. 
Since the Chern index of the |+) state is opposite to the 
ground state, this will introduce an error in the pumped 
charge proportional to the relative time spent in the ex- 
cited state. This is why high frequency noise must be 
thoroughly filtered and the pumped speed adjusted to 
quench Landau-Zener transitions. 

This study of quantization by controlled path, although 
conducted around CPP circuits, is quite general. The 
necessary criteria which have been stated in the intro- 
duction, are now discussed in more detail 

1. The adiabatic condition can only be satisfied for 
discrete spectra. Furthermore, adiabaticity only 
holds well for a quantum state if its energy split- 
tings with other levels are sufficient. 

2. The presence of isolated degeneracies between two 
lowest states states is required. For a complex 
Hamiltonian, this is possible only if it depends on 
three continuous parameters. In quantum circuits, 
tunable parameters are typically gate voltages and 
magnetic fluxes. For Cooper pair pumps, the pa- 
rameters are two gates charges and a phase (CPP) 



or two phases and a charge (cf. the sluice pump^^). 
The parameters are used to generate controlled 
paths in the three-dimensional parameter space P. 
Additional parameters are useful for adiabatic com- 
putations (section V). There are two recipes to lo- 
cate the degeneracies. The first one hinges on their 
topological signature on the Berry's phase (tt). As 
argued by A.J. Stone^^, one can divide the parame- 
ter space P in small volumes around which Berry's 
phase is computed. If a single degeneracy exists 
inside the loop, Berry's phase picks it up. Another 
method is to detect the vorticity of the eigenvector 
bundle in the vicinity of a degeneracy (see sections 
II and IV): in this case the "mirror symmetry" and 
"T-symmetry" are broken for the ground state, i.e. 
they flip the eigenvector bundles intersecting at the 
degeneracy (Kramers symmetry). 

3. Only observable proportional to the partial deriva- 
tive of the Hamiltonian H with respect to one of 
the tunable parameters i?^ {(p for the CPP) can be 
quantized using TQCP. In this case, its averaged 
expectation value over a closed surface iS around 
a degeneracy is proportional to the Chern index of 
this surface. When the variable is periodic, we may 
choose this surface with the topology of a torus. A 
path sweeping this surface or covering it densely is 
the controlled path expressing the quantization. 

Although topological quantization is quite robust, 
transitions to the first excited state (induced by fiuctu- 
ations or Landau-Zener processes) are problematic be- 
cause the Chern indices of the two lowest levels have op- 
posite sign. 

The above criteria only concern geometrical contribu- 
tions. But the dynamical evolution of the observable may 
also contribute. An accurate quantization is possible 
if they can be eliminated through symmetries or other 
schemes. This depends to some extend of the physical 
problem on hand. For quantum circuits, the following 
condition is sufficient to average them out: 

4. The Hamiltonian and the surface <S are periodic 
in the parameter ip, or have some symmetry with 
respect to if, such that the dynamical contribution 
averages to zero. 

The value of the topological quantization is its strong 
robustness to adiabatic parameters fluctuations, the 
Chern index of the surface enclosing the degeneracy be- 
ing the quantum number of the quantity of interest. For 
a Cooper pump, the magnitude of the current generated 

I = 2ecf\Ti)ue (61) 

is of order 30 pA for realistic values of the parameters. 

One of the potential application for a Cooper pair 
pump is as an accurate current source^° to close 
the metrology-triangle'^^, relating frequency to volt- 
age through the Josephson effect, voltage to current 
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through the quantum Hall effect and frequency to cur- 
rent through topological pumping. It is interesting to 
note that all these effects are the result of topological 
quantization^"^^. In order to be useful in this context, 
where the pumped current feeds a Hall bar, a higher cur- 
rent is needed (100 nA or higher''^). Large gain ampli- 
fication schemes are being designed to fulfill this condi- 
tion. Another method for a current-frequency conversion 
relies on Bloch oscillation'^''. This method has been re- 
cently demonstrated by the quantronium group'^^. This 
method, also based on a topological quantization (Thou- 
less pumping^^), also appears to be quite promising. 

To conclude, we have stressed the importance of topo- 
logical quantization in superconducting circuits. We 
showed how closed paths can be chosen in parameter 
space to generate a surface on which the "charge" Q is 



quantized by a Chcrn index. We have also shown how 
the charge transferred on any path could be obtained, 
giving also a local picture which cannot be derived from 
the usual global geometrical picture. Finally, the tech- 
nique developed around superconducting circuits is quite 
general and can be applied to any problems in which the 
criteria listed in the introduction and conclusion are sat- 
isfied. 
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